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Abstract. In this work we study the homogenization problem for nonlinear 
eigenvalues of quasilinear elliptic operators. We obtain an explicit order of 
convergence in k and in e for the (variational) eigenvalues. 



1. Introduction 

In this paper we study the asymptotic behavior (as e — > 0) of the eigenvalues of 
the following problems 



(1.1) 



-div(a e (x, Vu £ )) = X £ p e \u £ \ p ' 2 u £ in ft 
■ £ = on <9ft, 



where ft C is a bounded domain, e is a positive real number, and X 6 is the 
eigenvalue parameter. 

The weight functions p s (x) are assumed to be positive and uniformly bounded 
away from zero and infinity 

(1.2) < p~ < p £ (x) < p + < oo 

and the family of operators a e (x,^) have precise hypotheses that are stated below, 
but the prototypical example is 

(1.3) a E (x, Vu £ ) = A £ (x)\Vu £ \ p - 2 Vu £ , 

with 1 < p < +oo, and A £ {x) is a family of uniformly elliptic matrices (both in 
x e ft and in e > 0). 

The study of this type of problems have a long history due to its relevance in dif- 
ferent fields of applications. The problem of finding the asymptotic behavior of the 



eigenvalues of (1.1) is an important part of what is called Homogenization Theory. 
Homogenization Theory is applied in composite materials in which the physical pa- 
rameters such as conductivity and elasticity are oscillating. Homogenization Theory 
try to get a good approximation of the macroscopic behavior of the heterogeneous 
material by letting the parameter e — > 0. The main references for the homogeniza- 
tion theory of periodic structures are the books by Bensoussan-Lions-Papanicolaou 
|23j . Sanchez-Palencia [57], Olemik-Shamaev-Yosifian [23] among others. 



In the linear setting (i.e., a e (%,£) as in (1.3) with p = 2) this problem is well 



understood. It is known that, up to a subsequence, there exists a limit operator 



2010 Mathematics Subject Classification. 35B27, 35P15, 35P30. 

Key words and phrases. Eigenvalue homogenization, nonlinear eigenvalues, order of 
convergence. 



2 



J FERNANDEZ BONDER, J P PINASCO, A M SALORT 



cih{x, £) = A h (x)^ and a limit function p such that the spectrum of ( 1.1 ) converges 
to that of the limit problem. 



(1.4) 



—div(ah(x,S7u)) — Xp\u\ p 2 u in fi 
u = on <90, 



Let us recall that the convergence of eigenvalues in the multidimensional linear 
case was studied in 1976 by Boccardo and Marcellini [5] for general bounded matri- 
ces. Kesavan |21j studied the problem in an periodic frame. This is an important 
case of homogenization, i.e. when p £ (x) — p{x/e) and A e (x) — A(x/e) where p{x) 
and A(x) are Q— periodic functions, Q being the unit cube in M. N and in this case, 
the limit problem can be fully characterized and so the entire sequence e — > is 
convergent. See [2"Tll2"2"] . 

In the general nonlinear setting, Bafhco, Conca and Donato [4], relying on the 
G— convergence results of Chiado Piat, Dal Maso and Defranceschi [10 for mono- 



tone operators, study the convergence problem of the principal eigenvalue of ( 1.1 1. 
The concept of G— convergence of linear elliptic second order operators was intro- 
duced by Spagnolo in [29] ■ See Section 2 for the precise definitions. 

The convergence problem for higher (variational) nonlinear eigenvalues was ad- 
dressed by T. Champion and L. De Pascale in [§] where by means of the G— conver- 
gence methods they are able to show convergence of the (variational) eigenvalues 



of (1.1 1 to those of the limit equation (1.4 1 



Now, we turn our attention to the order of convergence of the eigenvalues that 
is the focus of this work. To this end, we restrict ourselves to the problems 



(1.5) 



I -div(a(x, W)Vu £ ) = X £ p e \u £ \P- 2 u e in fi 
\u £ =0 on dil, 

where the family of weight functions p e are given in terms of a single bounded 
Q— periodic function p in the form p s (x) := p(x/e), Q being the unit cube of M. . 

The limit problem is then given by 

J —div(a(x, Vm)Vu) = \p\u\ p ~ 2 u in fi 
I u = on dQ, 



(1.6) 



where p is the average of p in Q, i.e. P — Jq p{x) dx. 

The first estimate of the order of convergence of the eigenvalues, for the linear 
case, can be found in Chapter III, section 2 of [21] • By estimating the eigenvalues 
of the inverse operator, which is compact, and using tools from functional analysis 
in Hilbert spaces, it is prove that 

Here, C is a positive constant, and /J* satisfies 

< Pe < K\ 

with 

lim p k E = 

E-J-0 
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for each k > 1. 

The problem, again in the linear setting and in dimension N = 1, with a = 1, was 
recently studied by Castro and Zuazua in In those articles the authors, using 

the so-called WKB method which relays on asymptotic expansions of the solutions 
of the problem, and the explicit knowledge of the eigenfunctions and eigenvalues of 
the constant coefficient limit problem, proved 

\yM~V^\<Ck 2 e, 

equivalently, since the Weyl's formula implies that A| ~ k 2 , 

\X%-X k \<Ck 3 e 

and they also presented a variety of results on correctors for the eigenfunction 
approximation. Let us mention that their method needs higher regularity on the 
weight p, which must belong at least to C 2 (0) and that the bound holds for k ~ e -1 . 

In the linear problem, in any space dimension, Kenig, Lin and Shen [18 (allowing 
an e dcpendance in the diffusion matrix of the elliptic operator) proved that for 
Lipschitz domains £1 one has 

|A| - A fc | < Ce|log(e)|i+ CT 

for any a > 0, C depending on k and a. 

Moreover, the authors show that if the domain is more regular (C 1,1 is enough) 
they can get rid of the logarithmic term in the above estimate. However, no explicit 
dependance of C on k is obtained in that work. 

Later on, in |19j the authors obtain the precise dependence on k showing that 
|A| - A fc | < Cfc^e|log(e)|i+- 

for any a > 0, C depending on a. Again, when the domain Q is smooth, the 
logarithmic term can be removed. 



In this paper we analyze the order of convergence of eigenvalues of (1.5| to the 



ones of (1.6) and we prove that, 



|A% — A fc | <Ck^e 

with C independent of k and e. In this result, by A| and A& we refer to the 



variational eigenvalues of problems (1.5 1 and (1.6 1 respectively. 

Some remarks are in order: 

(1) Classical Weyl's estimates on the eigenvalues of second order A-dimensional 
problems show that X k and A| behaves like cfc« , with c depending only on 
the coefficients of the operator and N. Hence, the order of growth of the 
right-hand side in the estimate of [21] is 

A|(A fe ) 2 e5 kTf £ h a i 

— — — ~ < Ck N £ 2 . 

1 - \ h pi 1 - A fe /3| " 

Moreover, the constant involved in their bound is unknown. 
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(2) If we specialize our result to the one dimensional linear case, we recover 
the estimate obtained in [5]. Moreover, we are considering more general 
weights p since very low regularity is needed and the estimate is valid for 
any k. On the other hand, no corrector results are presented here. 

(3) In our result no regularity assumptions on fi are needed. Any bounded 
open set will do. 

(4) The constant C entering in our estimate, can be computed explicitly in 
terms of the weight p, the diffusion operator a(x, £) and the Poincare con- 



stant in the unit cube of K. (see Theorem 3.6 ) 



In the one dimensional problem, we can provide a extremely elementary proof 



of Theorems 3.5 and 3.6 generalizing the estimates obtained in [8]. Moreover, in 



this case an oscillating coefficient can be allowed. 

Let us mention that in the nonlinear case considered in this paper there are 
no quantitative estimates on the convergence of the eigenfunctions. However, it 
is possible to prove in one spatial dimension that the zeros of the eigenfunctions 
converge to the zeros of the ones of the limit problem and we also find an explicit 
estimate of the rate of convergence of the nodal domain. 



Organization of the paper. The rest of the paper is organized as follows: In 
Section 2, we collect some preliminary results on monotone operators that are 
needed in order to deal with (1.1) and some facts about the eigenvalue problem. 
We also recall the definition and some properties of G— convergence. In Section 
3 we prove the convergence of the fcth- variational eigenvalue of problem (1.5 1 to 
the fcth-variational eigenvalue of the limit problem (1.6), and we study the rates 
of convergence. In Section 4 we study the one dimensional problem, and we show 
in Section 5 that in this case, an oscillating coefficient can be allowed. In Section 6 
we deal with the zeros of eigenfunctions, and we close the paper in Section 7 with 
some numerical results. 



2. Preliminary results 

In this section we review some results gathered from the literature, enabling us 
to clearly state our results and making the paper self-contained. 

2.1. Monotone operators. We consider A: W ,p (Ct) — > W~ 1,p (Q) given by 

Au := — div(a(x, Vtt)), 

where a: fl x M. N — > R N satisfies, for every £ £ M. N and a.e. x € fi, the following 
conditions: 

(HO) measurability: a(-, •) is a Caratheodory function, i.e. a(x, •) is continuous 

a.e. igfi, and a(-,£) is measurable for every £ G R N . 
(HI) monotonicity: < (a(x, £i) — a(x, £2))(£i — £2)- 
(H2) coercivity: a|£| p < a(x, £)£. 
(H3) continuity: a(x,g) < /3|£| p_1 . 

(H4) p— homogeneity: a(x,t£) — i p ~ 1 a(cc,£) for every t > 0. 
(H5) oddness: a(x, — £) = — a(x, £). 
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Let us introduce ^(x, £1, £2) = a{x,£,i)£,i + a(x, £2)6 for all £1,^2 € R^, and all 
x £ ft; and let S = min{p/2, (p — 1)}. 

(H6) equi- continuity: 

\a(x, 6) - a(x, 6)| < c*(x, 6, 6) (p_1 - 4)/p (a(a:, 6) - a(x, - 6) 5/p 

(H7) cyclical monotonicity: Yli=i a {x, £,i)(£,i+i ~~ £«) < 0, for all k > 1, and 

6, • • • witn £1 = 6fe+l- 

(H8) strict monotonicity: let 7 = max(2,p), then 

See [4], Section 3.4 where a detailed discussion on the relation and implications 
of every condition (H0)-(H8) is given. 

In particular, under these conditions, we have the following Proposition: 

Proposition 2.1 ([3], Lemma 3.3). Given a(x,£) satisfying (H0)-(H8) there exists 
a unique Caratheodory function $ which is even, p— homogeneous strictly convex 
and differentiate in the variable £ satisfying 

(2.1) a\^<Hx^)<(3\C\ P 
for all £ G M. N a.e. x € ft such that 

V^(x,0=pa(x,0 
and normalized such that $(2;, 0) = 0. 

2.2. The nonlinear eigenvalue problem. In this section we review some prop- 
erties of the spectrum of for fixed e. That is, we study 

J — div(a(x, Vu)) = \p\u\ p ~ 2 u infl 
I u = on dft, 



(2.2) 



where < p_ < p(x) < p+ and < a < a{x) < (3 for some constants p—,p + ,a and 
P. 

We denote by 

S:={AeM: there exists it € Wj' p (f2), nontrivial solution to (2.2)} 

the spectrum of (2.2 1. 

By means of the critical point theory of Ljusternik-Schnirelmann it is straight- 
forward to obtain a discrete sequence of variational eigenvalues {Afe}fe £ N tending to 
+00. We denote by T, var sequence of variational eigenvalues. 

The fcth-variational eigenvalue is given by 

, La(x,\7v)\Vv\ 2 
2.3 A fc = inf sup Jo ; ' . . . ' 

cer kv £ J n p(x)\v\P 

where 

r fc = {C C W^ p {fl) : C compact, C = -C, 7(C) > k} 
and 7(C) is the Kranoselskii genus, see for the definition and properties of 7. 

For the one dimensional p— Laplace operator in J = (0,^), 
(2.4) - (\w'\ p - 2 w')' = p\w\ p - 2 w 



6 



J FERNANDEZ BONDER, J P PINASCO, A M SALORT 



with zero Dirichlet boundary conditions u>(0) = w(£) = 0, we have 



(2.5) 



inf sup 



Jj\v'\Pdx 



cec kve c Jj\v\Pdx' 

with v e Wq P (J). 

Here, all the eigenvalues and eigenf unctions can be found explicitly: 

Theorem 2.2 (Del Pino, Drabek and Manascvich, [H]). The eigenvalues /Jfc and 
eigenfunctions u>k of equation (2.4) on the interval J are given by 



Wk{x) = s\n p {i: p kx/t). 

Remark 2.3. It was proved in |15) that they coincide with the variational eigenvalues 
given by equation (2.5). However, let us observe that the notation is different in 
both papers. 

The function sin p (cc) is the solution of the initial value problem 

f-fU/lP-W = l7,|P- 2 7, 



and is defined implicitly as 



Moreover, its first zero is 7r p , given by 



v(0) = 0, ?/(0) = 1, 

(•sin* (a) / _ ^ 1 j p 



( — ) 



p- 1 \1 



1 - tP 



dt. 



dt. 



In problem (2.2) in one space dimension with a(x,u') = \u'\ p ~ 2 u' is studied 
and, among other things, it is proved that any eigenfunction associated to Afc has 
exactly k nodal domains. As a consequence of this fact follows the simplicity of 
every variational eigenvalue. 

The exact same proof of [3] works for genera a(x,u'), and so we obtain the 
following: 



Theorem 2.4. Every eigenfunction of (2.2) corresponding to the kth-eigenvalue 
\k has exactly k — 1 zeroes. Moreover, for every k, Afc is simple, consequently the 
eigenvalues are ordered as < Ai < A2 < • • • < Afc /*• +00. 

Now, using the same ideas as in [16] is easy to prove, for the one dimensional 
problem, that the spectrum of (2.2) coincides with the variational spectrum. In 
fact, we have: 



Theorem 2.5. S 



i.e., every solution of problem (2.4) is given by (2.5). 



Proof. The proof of this theorem is completely analogous to that of Theorem 1.1 

in pig. □ 
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2.3. Definition of G— convergence. We include the definition of G— convergence 
and the main results in [H [TU] for a sake of completeness, although we will need 
these facts only in Section 5. 

Definition 2.6. We say that the family of operators A £ u :— —div(a £ (x,'Vu)) 
G— converges to Au := — div(a(x, Vu)) if for every / <E W~~ l p (ft) and for every f £ 
strongly convergent to / in W~~ 1,p (51), the solutions u £ of the problem 

!—div(a £ (x,'Vu £ )) = f e in ft 
u £ = on dfl 

satisfy the following conditions 

u £ ->> u weakly in Wq' p (Q), 
a e (x, Vu £ ) a(x, Vu) weakly in (L p (ft)) N , 

where u is the solution to the equation 

{—div(a(x,Vu)) = f in ft 
u = on dfl. 

For instance, in the linear periodic case, the family — div(A(x/e)Vu) G-converges 
to a limit operator — div(A*'Vu) where A* is a constant matrix which can be char- 
acterized in terms of A and certain auxiliary functions. See for example [12] . 

It is shown in [3] that properties (H0)-(H8) are stable under G— convergence, i.e. 

Theorem 2.7 ( 4 , Theorem 2.3). If A £ u := —div(a £ (x,Vu)) G— converges to 
Au :— —div(a(x 1 Vu)) and a e (x,£) satisfies (H0)-(H8), then a(x,£) also satisfies 
(H0)-(H8). 

In the periodic case, i.e. when a £ (x,£) = a(x/e,£), and a(-,£) is Q— periodic for 
every £ £ 1 N , one has that A £ G— converges to the homogenized operator Ah given 
by AhU = — diti(o/,(Vu)), where a h : R N -> R N can be characterized by 

(2.6) a h (0= lim ~ [ a(z,Vxf + t)dx 

where £ £ K. , Q s (z s ) is the cube of side length s centered at z s for any family 
{z s } s >o in , and xt * s ^ ne solution of the following auxiliary problem 



(2.7) 



J -div(a(x,Vxt + t)) =0 in Q s {z s ) 
\xl £W^(Q s (z)), 

see [5] for the proof. 

In the general case, one has the following compactness result due to [TU] 

Proposition 2.8 ([10, Theorem 4.1). Assume that a £ (x,£) satisfies (H1)~(H3) 
then, up to a subsequence, A £ G— converges to a maximal monotone operator A 
whose coefficient a(x,^) also satisfies (H1)-(H3) 
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3. Rates of convergence 



In this section we prove that fcth-variational eigenvalue of problem (1.5) con- 
verges to the fcth-variational eigenvalue of the limit problem (1.6). 

Moreover, our goal is to estimate the rate of convergence between the eigenvalues. 
That is, we want to find explicit bounds for the error |A| — Afe|. 

We begin this section by proving some auxiliary results that are essential in the 
remaining of the paper. 

We first prove a couple of lemmas in order to prove Theorem |3.4| which is a 
generalization for p ^ 2 of Oleinik's Lemma |24j . 

We start with an easy Lemma that computes the Poincare constant on the cube 
of side e in terms of the Poincare constant of the unit cube. Although this result 
is well known and its proof follows directly by a change of variables, we choose to 
include it for the sake of completeness. 

Lemma 3.1. Let Q be the unit cube in R N and let c q be the Poincare constant in 
the unit cube in L q , q > 1. i.e. 

\\ u - ( u )q\\l«(Q) < c <jI|Vm|| L 9(q), for every u £ W 1,9 (Q), 

where {u)q is the average of u in Q. Then, for every u £ W ' q (Q £ ) we have 

\\u- (u) QE \\Li(Q e ) < c q e\\Vu\\ Lq{QE) , 

where Q £ = sQ. 

Proof. Let u £ W 1 ' q (Q e ). We can assume that (u)q c — 0. Now, if we denote 
u e (y) = u(ey), we have that u £ £ W 1,q (Q) and by the change of variables formula, 
we get 

|u|«= f \u e \ q e N <c q q e N ( \Vu e \ q = c q e q 

The proof is now complete. □ 
Remark 3.2. In pQ, Acosta and Duran show that for convex domains U, one has 

\\U~ (u)u\\ L i (U ) < |||Vu||ii((7), 

where d is the diameter of U. When we apply this result to the unit cube, we get 
(3.1) « < f . 

The next Lemma is the final ingredient in the estimate of Theorem |3.4| 



Lemma 3.3. Let O C l w be a bounded domain and denote by Q the unit cube in 
K. . Let g £ L oa (R N ) be a Q-periodic function such that 5 = 0. Then the inequality 



< llfllli«»(K^cie||Vu||ii(n) 



holds for every v £ W ' (O), where C\ is the Poincare constant given in (3.1). 
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Proof. Denote by I s the set of all z e Z N such that Q z , £ n ft ^ 0, Q z , e := e(z + Q). 



Given v G Wj' 1 (f2) we extended by outside fl and consider the function v s given 
by the formula 



v e (x) 



v(y)dy 



for x <E Qz.s- We denote by fii = U ze jeQ z e D f2. Then we have 



(3.2) 

Now, by Lema [3~T] we get 



9s(v - v E ) + 



\v - v e \\ L nn 



o = E 



Ida; 



(3.3) 



\Vv(x)\dx 



< cie||Vu||ii(n). 
Finally, since g — and since g is Q— periodic, we get 



(3.4) 



9eVe 



5> 



9e = 0. 



Qz.e 



Now, combining (3.3) and (3.4) we can bound (3.2) by 



9eV 



< 



l.9llL°°(R«)Ci£||Vw|| L i ( n). 



This finishes the proof. 



□ 



The next Theorem is essential to estimate the rate of convergence of the eigenval- 
ues since it allows us to replace an integral involving a rapidly oscillating function 
with one that involves its average in the unit cube. 

Theorem 3.4. Let g e L°°(1SL N ) be a Q— periodic function such that < g~ < g < 
g + < oo. Then, 

(9e{x) - g)\u\ p < pci\\g - 3|| i » 0(R w ) £||w||^ ( 1 n) ||Vu|| I/ » ) (n) 

for every u £ W ' P (Q). Here, c\ is the optimal constant in Poincare's inequality in 
L X (Q) given by (3.1). 



Proof. Let e > be fixed. Now, denote by h = g — g and so, by Lemma 3.3 
obtain 



(3.5) 

An easy computation shows that 

(3_6) IIV7n„,IPMI <r- vMn.WP- 1 



h e \u\ P < \\h\\L^)Cie\\V{\u\^)\\ L x m . 
lows that 

HvdOiUi^^PlK^IIVi/iu^o). 



Finally, combining (3.5) and (3.6) we obtain the desired result. 



□ 



Now we are ready to prove the main result of this section. 
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Theorem 3.5. Let A| be the kth-variational eigenvalue associated to equation (1.5) 
and let be X k be the kth-variational eigenvalue associated to the limit problem ( 1.6 I. 
Then 

|Af - A fe | < Cl p\\p-p\\oo{—) 1/p —e^{\ kl \l} 1+ -- 

a p 

Proof. Let 5 > and let G$ C Wq ,p (fi) be a compact, symmetric set of genus k 
such that 



A;. inf su P ^' Vu) - - 



sup J( \ ; _ ln ' +0(6). 



Ger kueG pJ n \u\P ueG * pJ n \u\P 

We use now the set G$, which is admissible in the variational characterization 
of the fcth-eigenvalue of ( |1.5| , in order to found a bound for it as follows, 



(3.7) 



X e k < sup 



/„<&(*, Vu) _ Vu) pJ Q \u\P 



4 IaPM p nee* P~Jn\ u \ P SnPM p 



*EGi Jo 



To bound A| we look for bounds of the two quotients in (J3/7J) . For every function 



u £ Gt we have that 



(3.8) 



2 < sup J " r y ' = X k + 0(5) 



Pin M p 



Since u € G k s C Wo' p (fi), by Theorem 



(3.9) 



Pfn H P 



< 1 + c lP e\\p- p\ 



Pfn\ v \ P 
re 

\\ u \\ P Lp(n)W Vu \\Lp(n) 



3.4 



we obtain that 

p-i 



Now, by ( |1.2[ ), (J2TJ) together with (|3_8|, we have 

"' f) , . 



P- V JqPM p 



(3.10) 



< 



p-V / n p|«|P y 

•pxl/p 1 



< {-)"" — (X k + 0(S))^ p . 
a p 



Then combining @, and ( |37l0| we find that 

AI < (A fe + 0(«5)) fl + c 1 p||p-p|| 00 (^) 1/p — e^ + O^)) 1 ^ 
V a P 



Letting 5 — > we get 
(3.11) 



A| - A fc < ClP ||p- p|U(£) 1/p ^sA^ + -. 



In a similar way, interchanging the roles of X k and A|, we obtain 
(3.12) 



A fe - A| < Cl p\\p-p\U^) 1/p -MXl) 1+1 ^ 

a p 



EIGENVALUE HOMOGENIZATION 



11 



So, from (3.11) and (3.12), we arrive at 



\K - \ k \ < c lP \\p~ p\\oo(— ) VP — £max{A fc , A £ fe } 1+ ?, 

a p 

and so the proof is complete. □ 

It would be desirable to give a precise rate of convergence in terms only con e, 
k and the coefficients of the problem. In order to achieve this goal, we need to 
give explicit bounds on the eigenvalues A& and A|. But this follows by comparison 
with the fcth-variational eigenvalue of the p— Laplacian, p^ and a refinement of the 
bound on p^ proved in |20j . 

In fact, from ( |2.1[ ) we have 

a f n \Vu\' < 4$(^Vu) < f Q \Vu\* 



P U U \ P ~ InP\ u \ P ~ P U U \ P ' 
from where it follows that 

max{Afc,A|} < — fi k . 

P 

Now, in [50], it is shown that 

_ / k -.p/N 
/**<Mi(p|) 

where /2i is the first Dirichlet eigenvalue for the p— laplacian in the unit cube. 
Finally, in [16] an estimate for jli by comparing with the first eigenvalue of the 
pseudo p— laplacian is obtained, namely 

£i < max{M p - 2 )/ 2 , 1}NttP, 

where ir„ is defined in Subsection 12.21 



Combining all of these facts, we immediately prove 

Theorem 3.6. With the same assumptions and notations as in Theorem \3.5\ we 
have 

|A fc -Af| <Cek P ~^, 
where the constant C is given by 

Remark 3.7. As we mentioned in the introduction, in the linear case and in one 
space dimension Castro and Zuazua [5] prove that, for k < Cs^ 1 , 

|A|-A fe | < Ck 3 e. 

If we specialize our result to this case, we get the same bound. The advantage of 
our method is that very low regularity on p is needed (only L°°). However, the 
method in 0, making use of the linearity of the problem, gives precise information 
about the behavior of the eigenfunctions u\. 
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Moreover, in the one dimensional linear case with diffusion coefficient equal to 
1, we can simplify the constant and obtain 

A fc - Afc < i — (nk) e 

{p r V p 

Remark 3.8. In [181 H9j- Kenig, Lin and Shen studied the linear case in any space 
dimension (allowing a periodic oscillation diffusion matrix) and prove the bound 

|A| - A & | < C* £ |lo ge | 1+CT A|. 

for some a > and C depending on a and k (The authors can get rid off the 
logarithmic term assuming more regularity on f2). If we specialize our result to this 
case, we cannot treat an e dependance on the operator, we get the same dependancc 
on k but without any regularity assumption on f2 we get the optimal dependence 
on e. 



(4.1) 



(4.2) 



4. The one dimensional problem 

In this section we consider the following one dimensional problems 

f-(|<r 2 <)' = A E p(f )krV in / := (0, 1) 
\u e {0) = u e (l) = 0, 

and the limit problem 

f-(KrM)' = A^(*)krv -mi 

\u e (0) = u e (l) = 0, 

where p is 1— periodic and p is the average of p in the unit interval. 

In order to prove the rate of convergence, let us assume first that p = and let 
R(x) = f£ p(t) dt. Then R is 1— periodic and if we denote R e (x) = R(^), we get 

sR' E (x)=p^). 

Hence, if v € W U (I), 

J p{~)v(x)dx = e [ R' e (x)v(x) dx = — e [ R e {x)v' (x) dx. 

So 



It is easy to see that 



i ,i 
p{f)v{x)dx <£||i?|U / \v'(x)\dx. 
o Jo 



|£|U<||max{p,0}||i = i||p||i, 



since p = 0. 

We have proved 

Lemma 4.1. Let v € W^ X {I) and p e L 1 ^) be such that p :— p{x) dx = 0. 
Then ^ 

r p^)v(x)dx < hpheWv'W,. 
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Then wc get the following corollary 
Corollary 4.2. Let u G Wq' p (I) and p G L X {I). The 



Proof. It follows just by noticing that 

(\u\ p y\ dx = P f 1 M'-V I < pIMirVll* 

Jo 

and applying Lemma |4.1| □ 



Now if we argue exactly as in Theorem |3.5| but use Corollary |4.2| instead of 
Theorem |3.4[ we get 

i 

|A|-A fc |<|||p-p||i^± £ max{A fc ,A|} 1+ ?. 



Ajfc,A| < — ^fc = — (7T p fe) P . 



The bound for Afc, A| follows directly from Theorem 2.2 In fact, 

1 1 

— A*fe = — 
P- P- 

So we have proved: 

Theorem 4.3. The following estimate holds 

|A|-A fc |<fJl^(^±)^)^. 

Remark 4.4. If we replace the unit interval I = (0, 1) by I( = (0,£) by a simple 
change of variables, the estimates of Theorem |4.3| are modified as 

(4.3) |A^(7,)-A fe (7,)|=F|A^(7)-A fe (/)|. 

5. The general equation 

In this section we consider, for the one dimensional problem the case where an 
oscillating coefficient in the equation is allowed, i.e., the problem 



~(o(f )W\ p - 2 u')' = A £ p(f )\u\p- 2 u in (0, 1) 
u (0) = u(l) = 



(5.1) 



We will show that this case can be reduced to Theorem |4.3| by a suitable change 
of variables. In fact, if we define 

and perform the change of variables 

(x,u) -> (y,v) 

where 

y = P e (x) = E P(f ), w(y)=u(a:). 
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By simple computations we get 

(-(\v\P-Hy = \tQ s (y)\v\P~ 2 v, y € [0,L e ] 
\v(0)=v(L e )=0 



where 
with 

and 



■ = d/dy, 
J„ o,(«)V(p-i 



ds — > L = a?- 1 . 



Qe(y) = a e ( x y/^p £ (x) 

= a{p-\l)fl^p{p-\l)) 

= Q(f)- 

Observe that Q is an L— periodic function. 

Moreover, it is easy to see that 
(5.2) \L e -L\<eL 
and that L e = L if e = 1/j for some j E N. 



In order to apply Theorem 4.3 we need to rescale to the unit interval. So we 
define 

w(z) = v(L £ z), z G / 

and get 

(-(\w\p- 2 w)- = LP\ £ Q e {L e z)\w\P- 2 w in / 
\w(0) = w(l) = 

So if we denote 5 = eL/L E , p s — LP\ £ and g(z) = Q(Lz) 7 we get that g is a 
1— periodic function and that w verifies 

{-{\w\p- 2 w)- = n s g(l)\w\ p ~ 2 w in / 
\w(0) = w(l) = 

Now we can apply Theorem 4.3 to the eigenvalues /i 5 to get 



(5.3) 



< 



P 115-9111 f9+\p 



2 a 1 - 



'6(ir p k) 



P+i 



In the case where e — 1/ j with j S N we directly obtain 

M-A fc |<^fi!^(^)Wr l . 

In the general case, one has to measure the defect between L and L e . So, 
(5-4) |A|-A fe | < l-Qti-Hkl + lim-Vl) < ^(U-Vk\ + j-K p p k p \LP-LP\). 



From (5.2) it is easy to see that 

\(^r~l\<p(l+s)P' 1 e. 
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SO 

(5.5) 



\L p - L p \ = L p \(^) p - 1| < P L p (l + e) p - 1 e. 



Finally, using ( |5.3[ ), ( |5.4[ ) and ( |5.5[ ) we obtain: 
Theorem 5.1. Let A| 6e t/ie kth— eigenvalue of 

f-Kf K)' = A £ p(f )| u r 2 « / = (o, i) 

|u(0) = u(l) = 

and Ze£ A& 6e i/ie kth— eigenvalue of the homogenized limit problem 

j-{a* p \u'\P- 2 u)' = Xp\u\P- 2 u ml 
1 u(0) = = 0. 



Then, if e = 1/j for some j € N, 

1 Pll5-3lli /5+ 



|A% — A*| < 



LP 2 gi V ff _ 



and if ' e =/= 1/j , 



|A€ — Afcl < 



i p ||g-g||i /g+\p g 



'e(7r p fc)P +1 



/3 



1 — £ p_ 



LP 2 «£ Kg 

6. Convergence of nodal domains 

In this section we prove the following result about the convergence of the nodal 
sets and of the zeros of the eigenf unctions. 



Theorem 6.1. Let (A|,w|) and (Afc,tifc) be eigenpairs associated to equations (4.1) 
and (4.2) respectively. We denote by Mf. and Nk to a nodal domains of u e k and 
respectively. Then 

-> |7Vfc| ase^O 

and we have the estimate 

\m\- p -Wk\- p \ < C£ (f +i + i) 



Proof. By using Theorem 3.6 together with (4.3) and the explicit form of the 
eigenvalues of the limit problem we obtain that 

(6.1) 

Also, 

(6.2) 



ai(/) = xum < Ai(A/fe) + cmrh < + C£ . 



k P TT P 

A|(/) > Afc(J) - cefc p+1 = - ce 



As Wk{x) = aia. v (kTr v x ) (se e Theorem 2.2) has fc nodal domain in / we must have 
|A4| =k~ 1 . Then by ([6j) and we get 



1 < 

cek p+1 < — ^ — p - + ce 



it follows that 
(6.3) 



Wk\- p - mr p < ce(k p+1 + 1) 
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Similarly we obtain that 



-p\U k \ 



Ai(A4) = A fe (J) > A|(J) - cek p+1 > A?(A/£) - cek p+1 



and using again Theorem 3.6 we get 

AfCA/fc) > Ai(A/"|) - ce 
it follows that 

(6.4) |A^|-f - |^4|"f < ce(kP +1 + 1). 



Ct 



Combining (6.3| and (6.4 1 the result follows. 



□ 



Finally, as a corollary of Theorem |6.1| we are able to prove the individual con- 
vergence of the zeroes of the eigenfunctions of (4.1) to those of the limit problem 
(p|. 



Corollary 6.2. Let (A|,u|) and (Xk,Uk) be eigenpairs associated to equations (4.1) 
and (4.2) respectively. Denote Xj and Xj, < j < k its respective zeroes. Then for 
each 1 < j < k 

Xj — > Xj when e — > 

and 

\xj-Xj] < ]ce{k p+1 + 1). 
In particular Xq = xq = and x\ = Xk = 1 by the boundary condition. 



Proof. With the notation of Theorem 6.1 we have that \Af%\ — > \Nk\- For the first 
pair of nodal domains we get 

\x\-xx\ =\xl-x e -x 1 +x \ = |K E a |-|A4 4 l| <ce(k p+1 + l) 

for the second couple 

\{xl - x 2 ) - (x\ ~ Xl )\ = \m, 2 \ - |A4, 2 || < ce(F +1 + 1) 

then 

\x £ 2 - x 2 \ < ce(k p+1 + 1) + \x\ -xi\< 2ce(k p+1 + 1). 



We iterate the reasoning for j < k, 
and the proof is complete 



x e j-Xj\ <jce(k p+1 + l) 



□ 



7. Some examples and numerical results 
We define the following Priifer transformation: 

, 71 \ \(l>=i) /P u ( x ) = P(x)S p {<p(x)), 

\u'(x) = p(x)C p (ip(x)) 
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As in f25j, we can see show that p(x) and <p(x) are continuously differentiable 
functions satisfying 

(72) jV(s) -(^)'+^|C P (^(x))r 2 C' P Mx))5 P M a; )) 
=j7$p{x)\S P (<p(x))\P 

and we obtain that 

Uk(x) = f "~~y^ pk(x)S p (ip k (x)), k>l 



is a eigenfunction of problem (2.4) corresponding to A& with zero Dirichlct boundary 
conditions. 



We propose the following algorithm to compute the eigenvalues of problem (2.4) 
based in the fact that the eigenfunction associate to A& has k nodal domain in I, 
so the phase function ip must vary between and kir p . It consists in a shooting 
method combined with a bisection algorithm (a Newton-Raphson version can be 
implemented too). 



Let a < A < b and let t be the tolerance 
Solve the ODE |7.2| and obtain tp\ and p\ 

Let w(x) = (p- I) 1 '? (\r(x))- 1/p p x (x)S p (vx(x)) 
Let a = w(l) 
while (\a\ > r) 
A = {a + b)/2 

Solve the ODE |7.2| and obtain tp\ and p\ 

Let w(x) = (p- l) 1 ^ (Xr(x)y 1/P p x (x)S p (^ x (x)) 
Let /3 = w(l) 
If {a(3 < 0) 

b= (a + b)/2 

else 

a= (a + b)/2 
end while 

Then A is the aproximation of eigenvalue with error < r 

For example, let us consider r(x) = 2 + sin(27ra;). In this case we obtain that 
r = J[2 + sm(2irx)dx — 2, and the eigenvalues of the limit problem are given by 

. l/p kiTp 

When e tends to zero the value of A e tends to the limit value A displaying oscilla- 
tions. 

When p = 2 the first limit eigenvalue is y% = tt/^/2 - 2.221441469. We see 
the oscillating behavior when plot y^Af as function of e in Figure 1 
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2,22145 

2 22 1 4 
222135 

2,2213 
2,22125 

2,2212 
2,22115 

2.2211 

0,002 0,004 0,006 0,008 0.01 0,012 0,014 0,010 0,018 0.02 

Figure 1 . The square root of the first eigenvalue as a function of 
e when r(x) — 2 + sin(27ra). 

A more complex behavior can be found in Figure [2] where we considered the 
weight r(x) = 2+si ^ 91TX - We observe that the sequence tends to 

Ai = tt 2 / / 1 dx = V3tt ~ 17.09465627. 

' J j 2 + sin 2ttx 

17.09B- 
17.095 
1-.35- ■ 

17.093 ■ ] 
17.09J- 

17.091 ■ \_ 

17.09 ■ 

17.089- " 

17.088 - 

17,087 -I 1 1 1 1 1 1 1 , 

0,1 0.2 0.3 0.4 0.E 0,8 0.7 0,8 

Figure 2 . The square root of the first eigenvalue as a function of 
e when r(x) = „ . - 1 „ — . 

It is not clear why the convergence of the first eigenvalue display the oscillations 
and the monotonicity observed (although the monotonicity is reversed for the weight 
r(x) = 2 — sin27ra;). We believe that some Sturmian type comparison theorem with 
integral inequalities for the weights (instead of point-wise inequalities as usual) is 
involved. However, we are not able to prove it, and for higher eigenvalues it is not 
clear what happens. 

Turning now to the eigenfunctions, with the weight r(x) — 2 + sin(27ra;), the 
normalized eigenfunction associated to the first eigenvalue of the limit problem is 
given by U\[x) = tt^ 1 sm(nx). Applying the numerical algorithm we obtain that 
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the graph of an eigenfunction associated to the first eigenvalue Af intertwine with 
the graph of U\(x). When e decreases, the number of crosses increases, and the 
amplitude of the difference between them decreases. In Figure [3] we can observe 
this behavior and the difference between u\ and u\ for different values of e. 




6=1 6=0.5 6=0.1 



Figure 3. The first eigenfunctions and the difference between 
them for different values of e. 

To our knowledge, it is not known any result about the number of the oscillations 
as e decreases, nor it is known if those oscillations disappear for e sufficiently small. 

The same behavior seems to hold for the higher eigenfunctions, see in Figure [4] 
the behavior of the fourth eigenfunction u\ when the parameter e decrease. 




e=l 6=0.5 e=0.1 



Figure 4. The fourth eigenfunctions and the difference between 
them for different values of e. 

Here, the convergence of the nodal domains and the fact that the restriction of an 
eigenfunction to one of its nodal domains Af coincides with the first eigenfunction 
of the problem in Af, together with the continuous dependence of the eigenfunctions 
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on the weight and the length of the domain, suggest that the presence or not of 
oscillations for the higher eigenfunctions must be the same as for the first one. 
However, the computations show very complex patterns in the oscillations. 
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